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We report on recent results for strongly frustrated quantum J1—J2 antiferromagnets in 
dimensionality d = 1,2,3 obtained by the coupled cluster method (CCM). We demon- 
strate that the CCM in high orders of approximation allows us to investigate quantum 
phase transitions driven by frustration and to discuss novel quantum ground states. In 
detail we consider the ground-state properties of (i) the Heisenberg spin- 1/2 antiferro- 
magnet on the cubic lattice in d = 1, 2, 3, and use the results for the energy, the sublattice 
magnetization and the spin stiffness as a benchmark test for the precision of the method; 
(ii) coupled frustrated spin chains (the quasi-one-dimensional J\— J2 model) and discuss 
the influence of the quantum fluctuations and the interchain coupling on the incommen- 
surate spiral state present in the classical model; (iii) the Shastry-Sutherland antifer- 
romagnet on the square lattice; and (iv) a stacked frustrated square-lattice Heisenberg 
antiferromagnet (the quasi-two-dimensional J1—J2 model), and discuss the influence of 
the interlayer coupling on the quantum paramagnetic ground-state phase that is present 
for the strictly two-dimensional model. 

Keywords: frustration, quantum antiferromagnets, quantum phase transitions, coupled 
cluster method 



1. Introduction 

Strongly frustrated quantum magnets have attracted much attention in recent years, 
both theoretically and experimentally. 1-3 In particular, quantum phase transitions 
have very much become the focus of interest (and see, e.g., Refs. [4-6]). At zero 
temperature, T = 0, there are no thermal fluctuations present and the transitions 
between ground-state (GS) phases are driven purely by the interplay of quantum- 
mechanical fluctuations and competition between interactions (e.g., frustration). In 
particular, novel quantum GS phases, such as valence-bond or spin-liquid phases 
may appear 6, which do not possess semiclassical magnetic long-range order (LRO). 
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A basic model which shows strong quantum fluctuations is the spin- 1/2 pure 
Heisenberg antiferromagnct (HAFM). The GS ordering of this quantum HAFM 
strongly depends on the dimensionality. 3 For example, unlike its three-dimensional 
(3d) and two-dimensional (2d) model counterparts on the cubic or square lattice, 
respectively, the one-dimensional (Id) pure HAFM (i.e., with nearest-neighbor an- 
tiferromagnetic bonds only, and all of equal strength) does not have a Neel-ordered 
GS. 6,8 Although the tendency to order is more pronounced in 3d quantum spin 
systems than in lower-dimensional ones, a magnetically disordered phase can also 
be observed for 2d and 3d HAFM's when strong frustration is present, e.g., for 
the J1-J2 model on the square lattice (see, e.g., Refs. [9-21] and references con- 
tained therein); or the 2d Shastry-Suthcrland antiferromagnct; 6,22-33 or, in three 
dimensions, for the HAFM on the pyrochlore lattice 34,35 or on the stacked kagome 
lattice. 36 

Besides the general interest in frustrated quantum antifcrromagncts as rich 
examples of quantum many-body systems, there is a strong motivation for their 
theoretical study which is driven by the many recent experiments on quasi-low- 
dimensional materials that are well described by a frustrated spin- 1/2 Heisenberg 
model. Among many others we mention here (i) the quasi-ld edge-sharing cop- 
per oxides like LICUVO4, LiCu 2 V02 and NaCu202, 37-40 which were identified as 
frustrated (with ferromagnetic nearest-neighbor (NN) and antiferromagnetic next- 
nearest-neighbor (NNN) interactions) quantum helimagnets with an incommensu- 
rate spiral GS; (ii) the quasi-2d SrCu 2 (1303)2, 41 ' 42 which is well described by the 
Shastry-Suthcrland antiferromagnet, and which exhibits a gapped quantum para- 
magnetic GS; and (hi) the quasi-2d Li 2 VOSi04 43,44 which can be described by the 
J\-Ji model on the square lattice. 

The theoretical treatment of the frustrated quantum antiferromagnets is far from 
being trivial. Though, surprisingly, one can find exact GS's of a simple product na- 
ture in some exceptional cases, 22 ' 45,46 many of the standard many-body methods, 
such as quantum Monte Carlo techniques, may fail or become computationally infea- 
sible to implement if frustration is present. Other methods, such as density- matrix 
renormalization group (DMRG) or exact diagonalization techniques, are essentially 
restricted to low-dimensional systems, at least for the present. Hence, there is con- 
siderable interest in any method that can deal with frustrated spin systems in 
any number of dimensions, including magnetic systems with incommensurate spiral 
GS's. A method fulfilling this requirement is the coupled cluster method (CCM). 
This approach, introduced many years ago by Coester and Kiimmel, 47 is one of the 
most universal and most powerful methods of quantum many-body theory (and for a 
review of which see, e.g., Ref. [48]). The CCM has previously been applied to various 
quantum spin systems with much success. 14 ' 33,49-62 The application to frustrated 
spin systems was started in the 1990's, 14,52 and has been developed in more recent 
years to the point where it has become a powerful tool in this field by including 
higher orders of approximations in a well-defined truncation scheme. 33 ' 54,56 ' 57,60 ' 62 

In this paper, we review some recent applications of the CCM to the generic 
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frustrated quantum spin-1/2 Hcisenberg model 

i,3 

where the Sj are the spin operators fulfilling sf =3/4, and the are the (compet- 
ing) exchange coupling parameters. 

The paper is organized as follows. In Sec. [2] we illustrate the main features of 
the CCM, paying special attention to the methodological aspects relevant for the 
application of the CCM to frustrated magnets. In Sec. [3] we report (as a benchmark 
test) high-order CCM results for standard unfrustrated lattices in d — 1,2,3 and 
compare them with other accurate methods. In Sec. [4] we discuss the influence of 
quantum fluctuations on GS spiral ordering for quasi-ld J\-Ji spin systems, and in 
Sec. [5]we report on a CCM treatment of the 2d Shastry-Sutherland model. Finally, 
Sec. [6] contains a discussion of the GS ordering of a stacked frustrated square-lattice 
HAFM (the quasi-2d J x -J 2 model). 

2. The Coupled Cluster Method 

The CCM formalism is now briefly outlined. For further details the interested reader 
is referred to Refs. [50,53-55,59-61]. The starting point for the CCM calculation 
is the choice of a normalized reference or model state |$), together with a set 
of (mutually commuting) multi-configurational creation operators {C^} and the 
corresponding set of their Hermitian adjoints {Cl}, 

<$|C+ =0 = C L |$) Vi^O, C+ = l; [C+,C+]=Q=[C L ,Cj]. (2) 

The operators C£ {Cl) arc defined over a complete set of many-body configurations 
denoted by the set of set-indices {L}. With the set {1$), } the CCM parametriza- 
tions of the exact ket and bra GS eigenvectors |\&) and (^1 of our many-body system 
are then given by 

\V)=e s \<P), S=^a L C+; (*| H$|£e- S , S=1 + $>lC l . (3) 

L#0 L^O 

The CCM correlation operators, S and S, contain the correlation coefficients, 
and cil, which have to be calculated. Once known, all GS properties of the many- 
body system can clearly be found in terms of them. To find the GS correlation 
coefficients and d^, we simply require that the GS energy expectation value 
H = (^iHl'fy) is a minimum with respect to the entire set {a/,, a/,}, which leads to 
the GS CCM ket-state and bra-state equations 

($|C^e- s ffe s |$) = , (^\Se- s [H,C+]e s \<i>) = VL ^ . (4) 

For frustrated spin systems an appropriate choice for the CCM model state |$) 
is often a classical spiral spin state, (i.e., pictorially, |$) = | \/*— *WS ■■■)), 
which is characterized by a pitch angle a. Such states include the Neel state, for 
which a = tt. In the quantum model the pitch angle may be modified by quantum 
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fluctuations. Hence, we do not choose the classical result for the pitch angle but, 
rather, we consider it as a free parameter in the CCM calculation, which is to be 
determined by minimization of the CCM GS energy. 

In order to find an appropriate set of creation operators it is convenient to 
perform a rotation of the local axes of each of the spins, such that all spins in the 
reference state align in the negative z-direction. This rotation by an appropriate 
angle Si of the spin on lattice site i is equivalent to the spin-operator transformation 

s* = cos <5isf + sin r^sf ; = ; s? = - sin f5,s? + cos f5jsf . (5) 

In this new set of local spin coordinates the reference state and the corresponding 
creation operators C~£ are given by 

I*) = I HU • • ■ > ; ct = j+ , sfsj , , . . . , (6) 

where the indices k, ... denote arbitrary lattice sites. In the new coordinates the 
initial Heisenberg Hamiltonian of Eq. (JTJ) becomes 

H = - J v { [cos(oty) + l](sfsj + §r§f) + [cos(aij) - l](s+s+ + §r§j) 

id 

+ 2sin(ay)[s 4 + s| - s-s+ + s~ s| - s?sj] + 4cos(ajj)s| s|} , (7) 

where ay = <5j — <5j is the angle between the two interacting spins, and sf = s^±is^. 
In the case of the Neel model state H becomes much simpler, so that, for example, 
for the case of NN bonds only, we have H = —JJ2(ij) (^t^f + l ^ r *7 ^~ 2sfs|). 

For the ensuing discussion of the GS properties we concentrate on the GS energy 
E, the order parameter m and the spin stiffness p s . Within the CCM scheme we 
have E = ($\e~ s He s \§) and m = -^(*| YliLi «f I*)- Tne s P in stiffness p s can be 
calculated by imposing a twist on the order parameter of a magnetically long-range 
ordered system along a given direction, i.e., 

E(9) E(0 = O) 1 

"~ /V - ~ ]v + 2^+°^)' W 

where E{0) is the GS energy as a function of the twist angle 9, and N is the number 
of sites (and where the interested reader is referred to Ref. [61] for further details). 

The CCM formalism is exact if we take into account all possible multispin config- 
urations in the correlation operators S and S. However, in general, this is impossible 
to do in practice for a quantum many-body system. Hence, it is necessary to use 
approximation schemes in order to truncate the expansions of S and S in Eq. ([3]) 
in any practical calculation. A quite general approximation scheme is the so-called 
SUBn-771 approximation. In this approximation all correlations in the correlation 
operators S and S are taken into account which span a range of no more than m 
contiguous sites and contain only n or fewer spins. In most cases the SUBn-n scheme 
is used (i.e., with n = to), and it is then (for spm-1/2 systems) called the LSUBn 
scheme. To find all different fundamental configurations entering S and S at a given 
level of LSUBrt approximation we use lattice symmetries and, where possible, any 
exact conservation laws. 
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Since the LSUBn approximation becomes exact in the limit n — > oo, it is useful 
to extrapolate the 'raw' LSUBn results to the limit n — ► oo. Although an exact 
scaling theory for the LSUBn results is not known, there is some empirical expe- 
rience 33,53-55 ' 62 indicating how the physical quantities for spin models might scale 
with n. For the GS energy we employ 54 ' 55 



Furthermore, we note that it may be useful to discard the LSUB2 results for the 
extrapolation, because generally they fit poorly to the asymptotic behavior. 55 For 
the order parameter and the stiffness one utilizes 55 ' 61 an extrapolation law with 
leading power 1/n, i.e., 



However, there is some experience that when applied to systems showing an order- 
disorder quantum phase transition this kind of extrapolation tends to overestimate 
the parameter region where magnetic LRO exists, i.e. to yield too large critical 
values for the exchange parameter driving the transition. 33 ' 54 ' 58 ' 61 The reason for 
such behavior might derive from the change of the scaling near a critical point. 
Hence, in addition to the extrapolation rule of Eq. (|10p for the order parameter m, 
we also use a leading 'power-law' extrapolation 33 ' 55,61 given by 



where the leading exponent C2 is determined directly from the LSUBn data. 

3. The CCM for the Pure HAFM on Cubic Lattices for d — 1, 2, 3 

During the last few years the running time and memory requirements of the original 
CCM code have been considerably improved. 60 ' 63 Consequently, it is now possible to 
run higher levels of approximation by using an improved parallelization procedure. 
In this Section we present a collection of CCM results for the (unfrustrated) pure 
HAFM (i.e., with the Hamiltonian of Eq. iQ} with J y = 1 for NN bonds, and J y = 
otherwise) on some basic lattices, and compare them with the most accurate results 
obtained by other methods. While some of the CCM results have already been 
published elsewhere, 55 ' 61 we also present here the new unpublished results from 
higher levels of approximation. The results are shown in Tables [TJ [2] and 02 in which 
the entries shown in boldface are the new ones. 

4. The Quasi-One-Dimensional J 1 J 2 Model 

One-dimensional quantum spin systems like the frustrated Ji-J2-model have at- 
tracted much attention over many years. 8 The physics of such quantum spin systems 






(10) 




(11) 
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Table 1. Data for the spin-1/2 linear-chain pure HAFM. Nf is the 
number of fundamental configurations for the Neel reference state, E/N 
is the GS energy per spin and m is the sublattice magnetization. The 
LSUBn results are extrapolated using Eq. (O for E/N and Eq. JToJl (in 
parentheses, Eq. JTTJ)for m. 



linear chain 


N f 


E/N 


m 


LSUB2 


1 


-0.41667 


0.33333 


LSUB4 


3 


-0.43627 


0.24839 


LSUB6 


9 


-0.44002 


0.20789 


LSUB8 


26 


-0.44137 


0.18297 


LSUB10 


81 


-0.44200 


0.16562 


LSUB12 


267 


-0.44234 


0.15263 


LSUB14 


931 


-0.44255 


0.14240 


LSUB16 


3362 


-0.44269 


0.13408 


Extrapolated CCM 




-0.44315 


0.07737 (-0.01086) 


Bethe ansatz bH 




-0.44315 


0.0 



Table 2. Data for the spin-1/2 square-lattice pure HAFM. Nf is the number of fundamental 
configurations for the Neel (in parenthesess, for the spiral) reference state, E/N is the GS 
energy per spin, m is the sublattice magnetization and p a is the spin stiffness. The LSUBn 
results are extrapolated using Eq. (J9j for E/N, Eq. j 10B . (in parentheses, Eq. llll ) for m and 
Eq. {TOjl for Pa- 



square lattice 


N S 


E/N 


m 


Pa 


LSUB2 


1 (3) 


-0.64833 


0.42071 


0.2574 


LSUB4 


7 (40) 


-0.66366 


0.38240 


0.2310 


LSUB6 


75 (828) 


-0.66700 


0.36364 


0.2176 


LSUB8 


1287 (21124) 


-0.66817 


0.35242 


0.2097 


LSUB10 


29605 (586787) 


-0.66870 


0.34483 




Extrapolated CCM 




-0.66936 


0.31024 (0.28073 ) 


0.1812 


3rd order SWT 64 




-0.66931 


0.3069 


0.1747 


QMC by 




-0.66944 


0.3070 


0.1852 



Table 3. Data for the spin-1/2 simple-cubic lattice pure HAFM. Nf is the number of 
fundamental configurations for the Neel (in parenthesess, for the spiral) reference state, E/N 
is the GS energy per spin, m is the sublattice magnetization and p a is the spin stiffness. The 
LSUBn results are extrapolated using Eq. (0 for E/N , Eq. I|10jl . (in parentheses, Eq. (TTJ) 
for m and Eq. (1 1 I t for p s . 



cubic lattice 


Nf 


E/N 


m 


Ps 


LSUB2 


1 (4) 


-0.89076 


0.45024 


0.2527 


LSUB4 


9 (106) 


-0.90043 


0.43392 


0.2416 


LSUB6 


181 (5706) 


-0.90180 


0.42860 


0.2380 


LSUB8 


8809 (444095) 


-0.90214 


0.42626 




Extrapolated CCM 




-0.90247 


0.42054 (0.42141) 


0.2312 


3rd order SWT bs 




-0.9025 


0.4227 


0.2343 



is often remarkably different from that of their corresponding classical counterparts, 
with a rich variety of different quantum GS's. As already mentioned in Sec.HJ novel 
edge-sharing copper oxides like LiCuV0 4 and NaCu 2 2 , 37-40 that were identified 
as quantum helimagnets with ferromagnetic NN and frustrating antiferromagnetic 
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T 



collinear state, (X =71 
h <0.25Ji 



Ji --Ji 



h 



spiral state, oc <7t 
h >0.25Ji 



Fig. 1. Illustration of the quasi-ld J1—J2 model and the reference states |4>) used for our CCM 
calculations. 



NNN bonds, have stimulated a great deal of renewed interest in frustrated Heisen- 
berg chains. In these materials a finite interchain coupling is also present, which 
in turn may lead to a low-temperature phase transition to a magnetically ordered 
phase. 

Therefore, we consider in this section the following quasi-ld spin-1/2 J\-J% 
Heiscnberg model, 

H = ( Jl X] Si '™ ' Si + 1 > n + ^ 2 X/ Si >™ ' S *+ 2 <™ J + Jj~y^ s i,n ' s i,n+l , (12) 

n \ i i / i,n 

where n labels the chains, J± is the interchain coupling, J\ is the in-chain NN 
coupling and J2 the in-chain NNN coupling, as shown in Fig. [TJ We fix J\ to J\ = 1 
(for the antifcrromagnetic case) or to J\ = — 1 (for the ferromagnetic case) and 
consider J 2 > (which is frustrating in both cases) and J± > 0. The case J± = 
corresponds to the strictly Id case. 

The classical system is characterized by a second-order transition from a collinear 
phase to a non-collinear spiral phase at |^-| = 0.25. For J2 > | J\ \ the classical spiral 
(pitch) angle a c \ is given by a c \ — arccos (— 0.25 J\j J2). Note that in the classical 
model neither the pitch angle a c \ nor the transition point \^-\ = 0.25 depends on 
the interchain coupling Ji . 

From the experimental point of view it is of interest to discuss the influence of 
quantum fluctuations on the pitch angle. Furthermore the question arises whether 
in the quantum model the interlayer coupling J± does or does not influence either 
or both of the pitch angle and the transition point. While the influence of quantum 
fluctuations on the pitch angle for J± = has been discussed previously in the 
literature, 52,69-71 this question has not been considered for finite J± > so far. 
Furthermore, any effect which a finite J± may have on the transition point has also 
not been considered in the literature. To discuss this point we use the CCM for the 
model of Eq. (fl2|) at the same level of approximation (viz., SUB2-3) as in Ref. [52]. 
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In that paper it was demonstrated that the SUB2-3 approximation for strictly Id 
systems leads to results of comparable accuracy to those obtained using the DMRG 
method. For the SUB2-3 approximation the relevant CCM equations can be found 
in closed analytical form, even for nonzero values of Jj_. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

J 2 J 2 



Fig. 2. The quantum pitch angle a qu versus J2 for the two cases (a) J\ = 1 and (b) J\ = —1 in 
the CCM SUB2-3 approximation, compared with the classical result. 




0.5 0.6 0.7 0.8 0.9 1 0.5 0.6 0.7 0.8 0.9 



a/71 a/71 

Fig. 3. Ground-state energy per spin in the CCM SUB2-3 approximation as a function of the 
pitch angle a for different (fixed) values of the NNN frustration strength J2 and antiferromagnetic 
NN bonds with Ji = 1, for the two cases (a) Jj_ = and (b) Jj_ = 0.7Ji. The quantum pitch 
angle a qu shown in Fig. [2] corresponds to the absolute minimum of the energy E. 

As quantum fluctuations may lead to a 'quantum' pitch angle a qu that is different 
from the classical value a c i, we consider the pitch angle in the reference state as 
a free parameter. We then determine a qu by minimizing the CCM energy E(a) 
with respect to a. Some results for the pitch angle a qu of the quantum model 
are presented in Fig. [2] We have checked explicitly that for J± = our results 
coincide with those of Ref. [52] . (We note that due to the different definitions of the 
Hamiltonian parameters that agreement is not obvious). For antiferromagnetic J\ 
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and weak interchain coupling Jj_ we find, contrary to the behavior of the classical 
system, a first-order transition from a collinear phase to a spiral phase, i.e., the 
quantum pitch angle jumps from a qu = 7r (in the collinear phase) to a qu < ir (in 
the spiral phase). This is illustrated in Fig. [3] by the energy curves E(a), which 
show a typical first-order scenario for J± = compared with a typical second-order 
scenario for J± = 0.7Ji. Furthermore, the transition point is both shifted to larger 
values of J2 > than in the classical case (and cf. Refs. [52,69,70]), and does now 
depend on Jj_ . For increasing Jj_ the transition point is shifted towards the classical 
value (but even for Jj_ = J\ it remains significantly above the classical value), 
and for large enough Jj_ the transition between the collinear and the spiral phase 
becomes continuous. 

For ferromagnetic J\ = — 1 one has, in complete analogy to the classical system, 
a second-order transition from the collinear to a spiral phase at J2 = 0.25 when 
J± = 0, which is in agreement with earlier considerations. 56 ' 71 However, for J2 > 
0.25| J\\ the quantum pitch angle deviates from the classical value and depends on 
the interlayer coupling J± . There is also a shift of the transition point for increasing 
values of J± towards larger values of J2. 

A common feature for both cases (i.e., J\ = — 1 and Ji = +1) is, that for 
increasing values of J2 the quantum pitch angle approaches its limiting value 7r/2 
much faster than for the classical model. 

The general behavior discussed above can be qualitatively related to the strength 
of quantum fluctuations (and see, e.g., the discussion in Rcfs. [33,54,56)], which 
themselves depend on the values of the exchange parameters of the model. In gen- 
eral, quantum fluctuations tend to stabilize collinear phases. Consequently, for the 
model with the strongest fluctuations, namely the strictly Id antiferromagnetic 
model (J± = 0,Ji = 1), the transition to the non-collinear GS takes place for 
largest J2 (cf. Fig. [2k). Increasing the interchain coupling then reduces the strength 
of the quantum fluctuations, which leads to a shift of the transition point to smaller 
values of J2. In the case of the strictly Id chain with ferromagnetic J\ = — 1 and 
J2 < 0.25|Ji|, the GS is the fully polarized ferromagnetic state, which does not 
exhibit any quantum fluctuations. As a result the transition point is the same as in 
the classical model, but the quantum pitch angle starts to deviate from the classical 
value immediately if J 2 > 0.25| Ji|. 

Finally, we emphasize that we did not speculate on possible spiral LRO. For 
strictly Id systems it is clear that there is no magnetic LRO and the quantum pitch 
angle discussed above corresponds to incommensurate short-range correlations. 

5. A Frustrated Heisenberg Antiferromagnet on the Square 
Lattice: The Shastry-Sutherland Model 

As already mentioned in Sec. [TJ the Shastry-Sutherland model, first introduced 
some 25 years ago, 22 has attracted much attention in connection with experiments 
on SrCu2(BOa)2. The model is characterized by a special arrangement of frustrating 
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NNN J2 bonds on the square lattice, as shown in Fig. 2] Its Hamiltonian reads 

h = Ji s, • Sj + j 2 ^2 s « ■ s fe ' ( 13 ) 

(id) {iM 

where the sum on runs over all NN bonds and the sum on {i, k} runs only over 
the selected NNN bonds shown in Fig.[U In what follows we set J\ — 1 and consider 




Fig. 4. Illustration of the Shastry-Sutherland model with antiferromagnetic NN bonds Ji = 1 
(solid lines) and NNN bonds J2 (dashed lines), together with its classical spiral state. The spin 
orientations at lattice sites n are given by the angles 8 = na c i, where n = 0, 1, 2, . . . , and a ci is 
the characteristic pitch angle of the classical spiral state. The state is shown for a c \ = 5n/6 and 
ra = 0,l,...,5. 



the case of (frustrating) J 2 > 0. Although the GS of this model is well understood 
in both the limits of small J2 and large J2, the GS phase at intermediate values of 
J2 is still a matter of discussion. 6,23-33 The CCM treatment of this model, briefly 
reviewed in this Section, is explained in more detail in Ref. [33]. 

We start with the classical GS of the Shastry-Sutherland model. It is the collinear 
Neel state for J2/J1 < 1, but a non-collinear spiral state for J2/J1 > 1 (and see 
Fig. 2] and Refs. [23,25]), with a characteristic pitch angle a c \ given by a c \ = ir for 
J% < J\ and a c \ = ir — arccos(Ji/ J2) for J2 > J±. The transition from the collinear 
Neel state to the non-collinear spiral state at J2I J\ = 1 is of second order for the 
classical model. We note further that there are only two different angles between 
interacting spins, namely a c i for the J\ couplings and — 2a c \ for the J2 couplings. 

Similarly as we did for the quasi-ld J\-J% model we now calculate the GS energy 
as a function of J 2 using as reference state a spiral state as shown in Fig. [4] Again, we 
find that due to quantum fluctuations the onset of the spiral phase in the quantum 
model is shifted to higher values of J2 , and the transition between the collinear and 
the spiral states becomes discontinuous, as seen from Fig. [5^. We conclude, that 
the collinear Neel state is the favored CCM reference state up to values J2 ~ 1.5Ji. 

It is known from the early paper of Shastry and Sutherland 22 that for large J2 
the quantum GS of the model of Eq. (flU)) is a rotationally-invariant orthogonal- 
dimer state given by |*)dimor = [I Ti)| ij) - I U)\ tj)]A/2 , where i and 
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Fig. 5. (a) The quantum pitch angle a qu for the Shastry-Sutherland model with J\ = 1 as a 
function of J2, calculated within the CCM-LSUBn approximation with n = 2,4,6; and (b) the 
GS energy per spin as a function of J2 of (i) the collinear Neel phase obtained within the CCM 
LSUBn approximation scheme with n = 4, 6, 8 and its extrapolated n — > 00 value using Eq. JS}, 
and (ii) the orthogonal-dimer state. 



j correspond to those sites which cover the J2 bonds. The energy per site of this 
state is E^ meI /N = — 3J2/8. We compare its energy in Fig. [SJa with that of the 
CCM GS obtained with the collinear Neel state as reference state. Our results 
demonstrate that the orthogonal-dimer state has lower energy than the Neel phase 
for J2 > 1.477 Ji, i.e., significantly before the point where the spiral CCM phase 
has lower energy than the collinear CCM phase. We note that we also checked that 
|^)dimor similarly remains the state of lowest energy in the region where the non- 
collincar spiral phase has lower energy than the Neel phase. We conclude that there 
is no spiral GS phase in the quantum model. The critical value J 2 = 1.477 Ji where 
the transition to the orthogonal-dimer phase takes place obtained by the CCM is 
in good agreement with results obtained by other methods (and see, e.g., Table 2 
in Ref. [31]). 

So far we have mainly discussed the energy of competing GS phases. The next 
question we would like to discuss is the question of the stability of the Neel LRO 
in the frustrated regime. For that purpose we calculate the order parameter (viz., 
the sublattice magnetization) m within the LSUBn approximation scheme up to 
n = 8 and extrapolate to n — > 00 using the extrapolation scheme of Eq. The 
results are shown in Fig. [6K- The extrapolated data clearly demonstrate that the 
LRO vanishes before the orthogonal-dimer state becomes the GS. The transition 
from Neel LRO to magnetic disorder is of second order. Hence, in agreement with 
previous investigations (and see, e.g., Ref. [31]), we come to the second important 
conclusion that there exists an intermediate magnetically disordered phase. The 
critical value Jf 1.14Ji where the Neel LRO breaks down agrees well with the 
corresponding value calculated by scries expansion techniques as given, for example, 
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in Table 2 of Ref. [31]. 
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Fig. 6. (a) The sublattice magnetization m versus J2 obtained within the CCM LSUBra ap- 
proximation with n = 4,6,8 and extrapolated to n — » 00 using Eq. lilt : and (b) the inverse 
susceptibility 1/xp versus J2 obtained within the CCM LSUBra approximation with n = 4,6,8 
and extrapolated to n — > 00 using Eq. 1101 . Inset: Pattern of plaquette valence bond order. 



The nature of the nonmagnetic GS in the region 1.14J2 < J2 < 1.477 J2 is a 
matter of some controversy in the literature. A favored candidate is a plaquette 
singlet phase. 25 ' 26 To address this question we follow similar reasoning to that used 
in Ref. [21] and consider the response of the spin system to a field F p given by 



x, y ' &w-\-l } y 



(14) 



p, namely the Hamiltonian of Eq. (|13p perturbed by the additional term of 

o21 



(15) 



where x,y are components (integer numbers) of the lattice vectors of the square 
lattice. This field corresponds to a plaquette valence bond order (as illustrated in 
the inset of Fig. [6b), which breaks the lattice symmetry. Thus, we use the CCM 
with the Neel state as reference state to calculate the energy per site e(Ji, J2, S) for 

Eq. (fT4|) . The susceptibility \p is then defined as 

_ ^1 

In Fig.[6jD we present the results for the inverse susceptibility, 1/xpj a s a function of 
J2. The extrapolation to n — > 00 is performed using an extrapolation scheme with 
leading power 1/n as in Eq. (flO|) . Clearly, in the magnetically ordered Neel phase 
X P is finite as it should be. However, close to the transition to the magnetically 
disordered phase at Jf « 1.14Ji the susceptibility becomes very large and diverges 
at J2 ~ 1.26Ji, which is close to Jf. Hence, we conclude from our CCM data 
that there exists a valence-bond phase between the Neel-ordered phase and the 
orthogonal-dimer phase. 
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6. The Quantum J\ J-2. Antiferromagnet on the Stacked Square 
Lattice 



As already mentioned in Sec. [T] the J1—J2 model on the square lattice is a canonical 
model to study quantum phase transitions in d — 2. 9-21 However, in experimental 
realizations of the J\-Ji model the magnetic couplings are expected to be not 
strictly 2d, since a nonzero interlayer coupling J± is always present. For example, 
recently Rosner et al. AA have found J± « 0.07Ji for I^VOSiCU, a material which 
can be described by a square lattice J1-J2 model with large J2. 43,44 

Therefore, in this Section we consider the influence of an interlayer coupling on 
the GS phases of the the J1-J2 spin-1/2 HAFM, i.e., we consider the HAFM on the 
stacked square lattice described by 



where n labels the layers and J± > is the interlayer coupling. The expression in 
parentheses represents the J\-Ji model of the layer n with intralayer NN bonds 
Ji = 1 and NNN bonds J 2 > 0. 

The classical GS's of the model are the Neel state for J2 < 0.5 Ji and another 
particular collinear state for Ji > 0.5Ji. The latter state (which we henceforth 
refer to as the collinear-columnar or, simply the collinear state) is a columnar (7r, 0) 
state characterized by a parallel spin orientation of nearest neighbors along the 
direction of one axis (say, the vertical or columnar direction) in each layer, and 
an antiparallel spin orientation of nearest neighbors along the perpendicular (say, 
horizontal or row) direction. It is well known 9-21 that for J± = the quantum model 
has two corresponding GS phases with semi-classical magnetic LRO, one (Neel -like) 
for small J2 < 0.4Ji and one (collinear-columnar- like) for large J2 > O.6J1, which 
are separated by a magnetically disordered (quantum paramagnetic) GS phase. 

For the treatment of the model of Eq. (fTB]) with arbitrary J± we apply the CCM 
and use both classical GS's (Neel and collinear-columnar) as reference states. Here 
we illustrate the CCM approach to this model only very briefly, and refer the in- 
terested reader to Ref. [62] for more details. In order to determine the GS phase 
transition points we calculate the order parameters for various values of J_l and de- 
termine those values J2 = 7Ncci(>/i) and J2 = jcoi(J±) where the order parameters 
vanish. In Fig. [7^ we present some typical curves showing the extrapolated order 
parameters (according to Eq. (fTT|) ) versus J2 for some values of J±. The magnetic 
order parameters of both magnetically long-range ordered phases vanish continu- 
ously as is typical for second-order transitions. We note, however, that there are 
arguments 11 that the transition from the collinear-columnar phase to the quantum 
paramagnetic phase should be of first order. The order parameters for both phases 
are monotonically increasing functions of J± , and the transition points 7Neei and 7coi 
also move together as J± increases. In Fig. [7b we present the dependence on J± of 
these transition points. Close to the strictly 2d case (i.e., for small J± <C Ji) the in- 




(16) 
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fluence of the interlayer coupling is largest. For a characteristic value of J*[_ m 0.19Ji 
the two transition points 7Ncei and 7 co i meet each other. Hence, we conclude that al- 
ready for quite weak interlayer coupling the magnetically disordered quantum phase 
that is present in the strictly 2d model disappears, and a direct transition between 
the two semi-classically ordered magnetic phases occurs. This conclusion is also in 
agreement with the statement that there is no magnetically disordered phase in the 
3d J\-Ji model on the bcc lattice. 72,73 




0.2 0.4 0.6 0.8 1 0.05 0.1 0.15 0.2 0.25 0.3 

Fig. 7. (a) The magnetic order parameter m versus Ji for various strengths of the interlayer 
coupling Jx (with Ji =1); and (b) the ground-state phase diagram (where the solid lines show 
those values of J2 for which the order parameters vanish). Note that Ji = 1. 
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